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Abstract 

Using finite basis sets, it is shown how to construct a local Hamiltonian, 
such that one of its infinitely many degenerate eigenf unctions is the ground- 
state full configuration interaction (FCI) wave function in that basis set. 
Formally, the local potential of this Hamiltonian is the optimized effective 
potential and the exact Kohn-Sham potential at the same time, because the 
FCI wave function yields the exact ground-state density and energy. It is 
not the aim of this paper to provide a new algorithm for obtaining FCI wave 
functions. A new potential is the goal. 
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Kohn Sham, Full Configuration Interaction 



1. Introduction 

In the Kohn-Sham model (KS) [H, 0], one considers a system of non- 
interacting fermions having the same ground state density, p(r), as a given 
system of interacting electrons. This is made possible by a convenient choice 
of the potential, u^s(r), in which the non-interacting fermions move. 

One can consider a model system of non-interacting fermions such that 
its ground state wave function minimizes the expectation value of the Hamil- 
tonian of a physical system of interacting electrons 0, 0]. The potential 
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describing this model system is the optimized effective potential, voEp( r ), 
and this also gives the name to the method (OEP). The OEP method has 
recently been applied to minimize energy expressions in the framework of 
KS. In this paper, however, we understand OEP exclusively as a method to 
minimize the expectation value of a physical Hamiltonian. 

With similar arguments like in References 0, @, 0] we construct a local 
potential for a finite orbital basis, such that all orbital energies are degener- 
ate. This allows us to choose the full configuration interaction (FCI) wave 
function as the ground state wave function. Thus, we obtain simultaneously 
the FCI density and energy. This means that the KS (exact density) and 
OEP (minimal energy) conditions can be satisfied within the given basis set. 

In this paper we point out that the FCI wave function is the ground state 
wave function of the local Hamiltonian constructed in this article if a finite 
basis set is employed. We are only interested in the potential. We do not 
suggest a new algorithm to obtain the FCI wave function. 

In section [2] we describe how to construct a fully degenerate non-local 
Hamiltonian. This is a non-local Hamiltonian with all orbitals degenerate. 
Furthermore, we show that in a finite basis a fully degenerate local Hamilto- 
nian can be constructed from the fully degenerate non-local Hamiltonian. In 
section [3] we analyze the conditions that must be satisfied to construct a fully 
degenerate local Hamiltonian. In section H] we give numerical examples for a 
fully degenerate local Hamiltonian. Section discusses the consequences of a 
fully degenerate local Hamiltonian. We argue that the FCI wave function is a 
ground state of the local Hamiltonian. In section [6] we draw the conclusions. 

2. Fully Degenerate Hamiltonians 

In this section we show, for a finite basis set, how to construct a local 
Hamiltonian that yields exactly the same orbitals and orbital energies as a 
given non-local Hamiltonian. The non-local Hamiltonian is fully degenerate, 
i.e. the orbitals are all degenerate. 

Assume a system with iV electrons and a finite basis set with M functions. 
We now search for a potential which makes all orbitals degenerate. Let 
us start with a local potential, v(r). This can be the nuclear potential, 
Vnuc{r) or this could be the KS or OEP potential, which include the nuclear 
potential. We obtain <fi m (r), the eigenfunctions of the non-interacting one- 
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particle Hamiltonian 

ho = - 1 -V 2 + v(r) (1) 
The eigenvalues are e m . We add a non-local potential 

M 



VNL 



^2 |0m) (e - e m ) (0m| (2) 



The eigenfunctions of the hamiltonian 

/ijvl = -^V 2 + v(r) (3) 

remain m (r), and they are all degenerate. Hence, it is possible to construct 
a Hamiltonian with a non-local potential that yields the same orbital energy 
for all orbitals. 

Furthermore, in a finite basis set, we can construct a local potential, 
vl(t), such that 

{<f>m \V L \ 4>n) = (4>m \v + V NL \ <f> n ) , Vm,B (4) 

Thus, 

h L = ~V 2 + v L (r) (5) 

produces the same (finite) hamiltonian matrix as h^L, and has all eigenvalues 
degenerate. Of course, the eigenstates of 

N 

H NL = J2 h NL^) (6) 

i 

and 

N 

H L = Y,h L (t) (7) 

i 

are also degenerate, since they are made of all the determinants constructed 
from the <p m . Thus, the ground state FCI wave function in the space spanned 
by the <p m is also an eigenstate of both H^l and Hl. As usual the eigenvalue 
of H L has no physical meaning. 
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3. Conditions for a Fully Degenerate Local Hamiltonian 

In the previous section we showed how to construct a fully degenerate 
local Hamiltonian, i.e. a local Hamiltonian with all orbitals degenerate. In 
this section we analyze the conditions that have to be met to insure that 
such a local Hamiltonian can be constructed. 



To obtain eigenvalues that are all degenerate it is necessary to satisfy all 
of the equations (jlj). In practice this can done by introducing a basis for the 
potential with, say K, functions. 



K 



v L (r) = v(r) + 2j hgt(r) (8) 
t 

For convenience we separated vl(t) in v(r), the local potential of equation ([TJ, 
and expanded only the non-local rest in {gt{r)}, the basis for the potential. 
Inserting this equation in the conditions (j4]) we obtain a set of equations 

K . 

J> / Mr)Mr)9t(r)dr = 5 kl (e - e k ) ,Vfc < I (9) 

for which solutions can be sought for b t . 

This is possible when equations ([9]) are consistent. However, linear de- 
pendencies in the products of orbitals can produce inconsistencies. Let us, 
thus, assume the linear dependence 

4>m{r)4>n{r) = ^2 c k i >mn (f) k (r)(f)i(r) (10) 

k<l 

We multiply each condition kl 7^ mn from ([9]) with the respective coefficient 
and sum them up. We obtain 



y^&t / ^ c k i >mn (j) k (r)(f)i(r)gt(r)dr = 8 kl ^ c k i, mn (e - e k ) (11) 

t J k<l k<l 
kl^mn kl^mn 

We have now used all condtions of equations (jSJ) but the one for the pair mn. 
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We subtract the remaining condition to obtain 



A' 



t 

( 



\ 



>ki 



k<l 
\kl^mn 



9t{r)dr 



E 

k<l 
\kl^mr 



CkL; 



;i2) 



Now we use the assumed linear dependency (cf. ([TO]) ) to obtain 



= 5 ki 



ki 



i 

E 

k<l 
\kl^mr 



\ 



CkL' 



(e - e k 



(13) 



/ 



If we assume for a moment, that none of the products (pi4>i are involved 
in the linear dependency of equation (ITUj) . the final equation ( TT3"j) is trivially 
satisfied. In this case m ^ n and Cki t mn = 0, Vfc = /, and hence the right 
hand side is equal to zero. 

If only one product is linearly dependent we obtain = Si — e. 
This equation can be satisfied by choosing e = £j. However, if more than 
one product 0j</>j is linearly dependent on the rest of all orbital products, 
equation f[T3l can not be satisfied. 

In this section we found that the diagonal terms (fiicfii must be linearly 
independent from the rest of all orbital products to guarantee a fully degen- 
erate local Hamiltonian. 



4. Numerical Example 

In this section we present numerical examples for a fully degenerate, finite 
basis, local Hamiltonian, i.e. a local Hamiltonian with all orbitals (finite 
number) degenerate. We consider He together with cc-pVDZ, cc-pVTZ and 
cc-pVQZ as an orbital basis, {Xi}, where we removed d- and f-functions. For 
the potential basis, {gt}, we choose an even tempered basis of s-type gaussian 
orbitals with exponents ranging from 0.1 to 40 000. For the cc-pVDZ orbital 
basis we use 5 potential basis functions, for the cc-pVTZ we use 10 potential 
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basis functions and for cc-pVQZ we use 25 potential basis functions. We use 
more potential basis functions than needed to insure that the space spanned 
by the orbital products is covered by our potential basis. The singular value 
decompositon used in our implementation will pick out only the necessary 
potential functions. 

The goal is to construct a Hamiltonian matrix, which has all eigenvalues 
equal, i.e., there is an orbital basis, in which it is equal to some e times the 
unit matrix. We choose e — 1 for convenience. In a first step we calculate 
the one-electron matrix, T, and the potential matrices, Gt, in the orbital 
basis. 

1 2 2 
--V 2 - - 

2 r 



T H = ( X 



Xi (14) 



Gt,ij = (xi\gt\xj) (15) 

The goal is to find a set of b t , such that 



K 



I = T + J2b t -Gt (16) 

t 

holds, where I is the unit matrix. Collecting the matrices Gt in a super 
matrix G, where the G t form the t-th column and rearranging equation (Fl6l) 
yields clS £L solution 

b = (I-T)-G- 1 (17) 

Since the matrix G is, in general, not square but rectangular, we use a 
singular value decomposition to calculate the pseudo- inverse G -1 . 

With the cc-pVDZ orbital basis we are able to obtain the unit matrix 
with an accuracy of 10 -14 . For the cc-pVTZ orbital basis the accuracy is 
1CT 11 and for cc-pVQZ the accuracy is 10 -7 . Figures [Q [2] and [3] display the 
expanded part of the potentials (cf. equation (jSj)) that yield the unit matrix 
in the given orbital basis. Strong oscillations are found close to the nucleus, 
as was reported in similar calculations @] . In comparison we also show a 
very accurate KS potential in Figure @] [8| . 

When comparing the corrections to the electron-nucleus potential we no- 
tice that vl having the FCI wave function as a solution (figures [1] - [3]) is 
quite different from the accurate KS one (figure @J. Not only the shape 
largely differs, but also the order of magnitude (Please notice the change of 
scale between figures.). Furthermore, no convergence towards the correct KS 
potential can be seen. 
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5. Discussion 



In this paper we want to obtain the OEP and KS potential for He in a 
given finite basis set. The KS potential is that potential that yields the FCI 
density 0] • In this respect the FCI density must be known to determine the 
KS potential. The OEP potential, on the other hand, is that potential whose 
ground state wave function minimizes the expectation value of the physical 
Hamiltonian 

In general, vks{ t ) and voEp{f) are not identical. It is well known, that 
the OEP potential of Helium, for a complete basis, is the Hartree potential. 
At the same time it is clear that the KS potential must differ, since the 
Hartree potential does not yield the exact density. 

Both, for the KS and the OEP models, there is no interaction between 
fermions and thus, in general, the ground state of the system can be described 
by a single Slater determinant. In the case of degeneracy of two or more single 
Slater determinants, any linear combination of the degenerate single Slater 
determinants is also a ground state. 

In the previous sections we showed how to construct a fully degenerate 
local Hamiltonian for a finite basis set. As a consequence, each and every 
wave function is a ground state of the local Hamiltonian. The OEP procedure 
demands to pick that wave function that minimizes the expectation value of 
the physical Hamiltonian. Doubtless, this must be the FCI wave function. 

There is no wave function that yields a lower energy than the FCI wave 
function. Consequently, the potential constructed in the previous sections 
is the OEP potential in the given finite basis set. At the same time this 
potential yields the FCI density. Consequently, the potential is also the KS 
potential. 

The potential that we constructed is not unique. It will differ if a different 
potential basis gt{r) will be used. The choice of s, the energy at which all 
orbitals are degenerate, will also influence the potential. Finally, the exact 
density can also be obtained from the local potential of a Schrodinger-like 



equation for yj p{r) [1 

To obtain the OEP or KS wave function it does not suffice to solve equa- 
tion (fT7|) and diagonalize the corresponding local Hamiltonian. More effort is 



needed. The FCI wave function must be construced in the usual way [12 
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6. Conclusion 



In the optimized effective potential (OEP) and the Kohn-Sham (KS) 
models one usually tries to avoid degeneracies. In this paper, however, we 
focus on the consequences of degeneracies. In an approach similar to the ones 
taken in References [1, @] we construct a fully degenerate local Hamiltonian. 
This means that all orbitals have the same energy. 

The independent particle Hamiltonian constructed with Vl(t) (cf. eq. 
(j3J)) simultaneously yields the ground state density for this basis set (the FCI 
density) and the lowest possible expectation value of the physical, interacting 
Hamiltonian (the FCI energy). Thus, the model Hamiltonian corresponds to 
both the KS and OEP solutions. Consequently, we have construced a KS 
and OEP potential at the same time. We did not construct the FCI wave 



function, which may be obtained from standard procedures [ll|, ll2 



We believe that the conditions outlined in section [3] can not be satisfied 
with very large basis sets. In section H] we see already for the cc-pVQZ basis a 
deviation of 10™ 7 . We conjecture, that equations (JHJ) cannot be satisfied, when 
the size of the basis set increases, due to the over-completness of the orbital 
products basis. This has also been noticed when attempting to construct the 
density matrix from the density 0, 0] • 
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Figure 1: The displayed potential is the local potential minus the nuclear potential, vl{t) — 
Vnuc{r) (for the potential basis see text). It yields the unit matrix in the cc-pvdz orbital 
basis. 
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Figure 2: The displayed potential is the local potential minus the nuclear potential, — 
Vnuc{r) (for the potential basis see text). It yields the unit matrix in the cc-pvtz orbital 
basis. 
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Figure 3: The displayed potential is the local potential minus the nuclear potential, vl{t) — 
Vnuc{r) (for the potential basis see text). It yields the unit matrix in the cc-pvqz orbital 
basis. 



12 




Figure 4: A very accurate KS potential for He minus the nuclear potential. 
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